Atomic matter-wave revivals with definite atom number in an optical lattice 
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We study the collapse and revival of interference patterns in the momentum distribution of atoms 
in optical lattices, using a projection technique to properly account for the fixed total number of 
atoms in the system. We consider the common experimental situation in which weakly interacting 
bosons are loaded into a shallow lattice, which is suddenly made deep. The collapse and revival 
of peaks in the momentum distribution is then driven by interactions in a lattice with essentially 
no tunnelling. The projection technique allows to us to treat inhomogeneous (trapped) systems 
exactly in the case that non-interacting bosons are loaded into the system initially, and we use 
time-dependent density matrix renormalization group techniques to study the system in the case of 
finite tunnelling in the lattice and finite initial interactions. For systems of more than a few sites and 
particles, we find good agreement with results calculated via a naive approach, in which the state 
at each lattice site is described by a coherent state in the particle occupation number. However, 
for systems on the order of 10 lattice sites, we find experimentally measurable discrepancies to the 
results predicted by this standard approach. 



I. INTRODUCTION 



In the early days of Bose Einstein condensates (BECs) 
with atomic gases there was a hvely debate how to recon- 
cile the notion of the phase of a BEC with atom number 
conservation [li3 . In condensed matter physics this dis- 
cussion goes back to Leggett's work of spontaneously bro- 
ken gauge symmetries, and the observabiUty and mean- 
ing of an absolute phase of a condensate, and of a relative 
phase between two BECs [I] [5] . An absohite condensate 
phase - provided we understand the BEC order param- 
eter t)) as expectation value of the boson particle 
destruction operator ~ implies a coherent superposition 
of states of atom numbers as in a coherent state. How- 
ever, in a theory where all observables commute with the 
number operator, the phase of such superpositions re- 
mains unobservable on a fundamental level. On the other 
hand, the relative phase of two condensates, associated 
with superposition states of atom number differences be- 
tween two BECs, is consistent with particle number con- 
servation, and is observable in an interference experiment 
which measures a correlation function (■0'''(x, t)Vj(x', t')) 
[3] . As in the case of coherence of light in quantum optics 
[S] , the factorization property of the correlation function 
(V''''(x, t)^(x', i')) — >• tp* {x,t)tp{-x.' ,t') for a given quantum 
state implies full visibility of the interference pattern, and 
provides a measurement of the BEC order parameter in 
a number conserving context. The emergence of inter- 
ference fringes in a quantum measurement process based 
on successive detection of atoms for two BECs prepared 
initially in a product of Fock states (eigenstates of fixed 
particle number) was analyzed in a language and formal- 
ism familiar from photon detection in quantum optics in 
Rcfs. [3 . Detection of an atom in an interference exper- 
iment erases the "which path" information from which 
BEC the atom came from, and thus prepares a superpo- 



sition state of atom number differences compatible with 
number conservation. 

Collapse and revival dynamics due to interactions in 
a BEC loaded into an optical lattice, as described by a 
Hubbard model and observed in the visibility of inter- 
ference fringes after releasing atoms from the trap [7] , 
is a second example where traditionally atomic number 
superpositions in the form of coherent states are invoked 
in the theoretical derivation and interpretation [SUllj. 
Somewhat surprisingly, a derivation of this phenomenon 
in a discussion respecting number conservation seems to 
be missing in the literature. Below we will provide an- 
alytical collapse and revival formulas for finite particle 
number, where we see the emergence of the usual re- 
sults in the limit of an infinite number of wells. On the 
more practical side, we will also connect the collapse and 
revival signal with a number conserving time-dependent 
DMRG treatment of Bose- Hubbard dynamics in a ID op- 
tical lattice. This makes it possible to calculate system- 
atically corrections due to finite tunneling between the 
wells and finite initial interactions which has not been 
possible in previous work. 

Collapse and revival of interference patterns driven by 
two-body interactions was first raised in discussing the 
dynamics of Bose-Einstein Condensates (BECs) in a dou- 
ble well [ini E], and later was investigated with cold 
atoms in optical lattices [S [7]. In these lattice exper- 
iments, weakly interacting bosons are typically loaded 
into a shallow lattice, which is then suddenly made deep. 
This prevents tunnelling of atoms between neighbouring 
sites, and the interference pattern corresponding to the 
momentum or quasi-momentum distribution of atoms in 
the lattice undergoes collapses and revivals due to on-site 
interactions. Most recently, these dynamics were used to 
investigate the dependence of on-site interaction energy 
shifts on the occupation number of a lattice site in a 3D 
optical lattice [3 [T^] . Dynamics of collapses and revivals 
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driven by two-body interactions have also been studied in 
the context of hard-core bosons superlattice [13], Bloch 
oscillations in a tilted lattice [TU [TS], and double well 
superlattice experiments, ]1M 117) . 

In essentially all of these cases, the experiments have 
been described based on an ansatz in which the initial 
state of the non-interacting BEC can naively be repre- 
sented as a product of local coherent states in the onsite 
occupation number. This can be represented as 



IBEC, 



M 

n- 



vac) 



(1) 



The rest of this paper is organized as follows: In Sec, [ill 
we first present the details of the system we study, includ- 
ing the Bose-Hubbard model, and the quasi-momentum 
distributions and correlation functions we calculate. In 



Sec. Ill we present exact results for revivals in this sys- 
tem, accounting properly for the conserved total number 
of particles in the system, and comparing the results to 
the naive approach outlined above. We account also for 
trapping potentials, finite tunneling in the lattice, and 
interactions in the initial state. In Sec. |IV| we present a 
summary and outlook for this work. 



where, denotes the creation operator for particles on 
site i, and is the local mean particle number. As 
we show below, this ansatz reduces to the correct non- 
interacting ground state in an infinite system, but in a 
finite-size system does not have a fixed total number of 
particles. 

In a deep lattice, the time-evolution of this state is 
simply governed by two-body onsite interactions gener- 
ated by the Hamiltonian term {U/2) niini — 1), with 
the number operator fii = b\bi. We find that expectation 
values such as the single-particle density matrix (SPDM) 
factorize, so that for a homogeneous system, we obtain 

(;.^i): 



{b\b,)t^{h);{b,)t 



= |/3|2gl/3r(2cos(;7t)-2) 



(2) 
(3) 



with Pi = (3 m. the homogeneous system. 

This naive number non-conserving ansatz leads to 
wrong results for small particle numbers, which could 
be quantitatively measured in an experiment. In this pa- 
per, we demonstrate this by calculating the exact results 
for finite systems and a fixed total number of particles 
via a projection operator approach. At the same time, 
we confirm that the naive ansatz yields accurate results, 
even for suprisingly moderate-sized systems larger than 
about five particles on five lattice sites. 

The approach we use here is similar to that used in pre- 
vious studies of Bose-Einstein condensates, where a Fock 
state of fixed particle number is written as a phase aver- 
aged coherent state [3] . Our results for a finite system are 
in agreement with previous work for homogeneous sys- 
tems based on a binomial expansion (see Ref. [E]), and 
we are also able to treat systems with open boundary con- 
ditions and trapping potentials. Using an optimised ver- 
sion of the Time-Evolving Block Decimation algorithm 
(TEBD) [m [20], we also study and discuss the effects 
of finite tunnelling in the lattice during time evolution 
pn 122) . as well as weak interactions in the initial state. 
Though for clarity we focus on results here in which the 
atoms are arranged along one dimension (e.g., with tight 
confinement in a 3D optical lattice), the results can be 
generalized straightforwardly to higher dimensions. 



II. THE MODEL 

We consider cold atoms loaded into the lowest Bloch 
band of an optical lattice , which are described by 

the Bose-Hubbard Hamiltonian {h = 1) 

U[t) 

(4) 



H 



-j{t) T b!b 



where are the bosonic annihilation (creation) opera- 
tors at site i, J denotes the tunneling amplitude between 
neighboring sites, U denotes the onsite interaction energy 
shift for two atoms, fii = blbi and Si is the energy offset 
on different sites, which arises due to external trapping. 
We will assume that we have a total of N particles on 
the M lattice sites. 

We assume that at time t = 0, the system is loaded in 
the ground state of the Hamiltonian in a relatively shal- 
low lattice where that atoms are weakly interacting. At 
time t — the depth of the lattice is rapidly increased, 
on a timescale that is fast compared with the tunnel- 
ing time ^ but slow compared with excitations to 
higher Bloch bands. In this way, the parameters J and 
U are changed suddenly at time t = 0, so that 



Jit) 



Jo, t<0; 
J t>0. 



U(t) 



Ua, t<0; 
U t>0. 



(5) 



Below we will first analyze the ideal case, of an initial 
non-interacting BEC, so that Uq = 0, and we have N 
particles all in the same single-particle ground state of the 
Bose-Hubbard model. For a single particle, the ground 
state in this system can be written as 
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M 

E 



/3i&^|vac), 



(6) 



with complex amplitudes /3i that are dependent on the 
external potential e^. Note that we normalize these am- 
plitude to the particle number A^, i.e. = A^, for 
reasons which become clear below. Then, the initial state 
for A^ particles reads 



IBECjv) = 



1 



f M ^ 



vac). 



(7) 



3 



In the ideal case, after increase of the lattice depth, the 
tunneling will be negligible on the timescales of interest 
after the lattice depth is increased, i.e., J « 0. Both of 
these ideal conditions will be relaxed in Sec. IIII El 

As mentioned above, collapses and revivals as studied 
in the experiments are observed via modulation of the 
visibility of interference fringes in the density distribu- 
tion of the atomic cloud after a free expansion (time-of- 
flight measurement) . This effectively measures the quasi- 
momentum distribution Uq {t) after a certain hold time in 
the lattice t, defined as: 



1 



M 



(8) 



where the discrete site-positions in the lattice, and 

{b\bj)t is the time-dependent single particle density ma- 
trix (SPDM). Note that so far we consider a system in 
arbitrary dimensions and with arbitrary lattice geome- 
try. Thus, in general the Xi and q are vectors. Along a 
single dimension, we consider a lattice constant a, and 
the discrete positions in that dimension can be written 
as x}^ = ai. Thus, along that dimension, which consists 
of AfiD sites, the quasi-momenta can have the discrete 
values q^^ = n2n/aM^^ (with integer n = l...M^^). 
In the following we will focus on the calculation of the 
time-evolution of the SPDM elements. 



III. EXACT RESULTS FOR REVIVALS IN THE 
BOSE-HUBBARD MODEL 

In this section, we present our results for revivals in the 
Bose-Hubbard model, correctly accounting for the fixed 
total number of particles in the system. We begin with 
a simple example of where the naive treatment breaks 
down, by studying two particles on two lattice sites in 
Sec. Ill A[ In Sec. |III B we then provide technical details 
of the projection method used to derive exact results in 
the case that atoms are initially non-interacting {Uq — 0) 
and that tunneling in the deep lattice is negligible (J = 
0). The exact results for these cases are then summarized 
and compared with the naive approach for homogeneous 
systems in Sec. |III C and for inhomogeneous systems in 
Sec. Imp] In Sec. |IIIE| we will relax these conditions 



(Uq — J = 0), accounting for finite timnelling and finite 
interactions in the initial state by performing essentially 
exact many-body numerical computations. 



A. Two particles on two lattice sites 

We begin with a simple example for which the ansatz 
of on-site coherent states, Eq. ([!]) clearly gives the wrong 
results. This can be easily demonstrated for two parti- 
cles on two lattice sites {N = M = 2) without exter- 
nal trapping, for which the initial state will be given by 
\^2)o = (|02) -I- |20))/2 + \n)/V2, where \n1n2) denotes 



a state with ni particles on the first site and 712 particles 
on the second site. Time-evolving this state under the 
Hamiltonian Q gives 

IV'2)t = I (e-**|02) + e-^*|20» + (9) 

so that the off-diagonal SPDM elements in a two-site 
system with two particles are given by 



{b\b2) = (blbi) ^ cos (Ut). 



(10) 



This clearly contradicts the result from Eq. as can 
be seen in Fig. [T] 



B. Particle number projection method 

In order to conveniently calculate the time evolution 
for an arbitrary number of particles, we now introduce an 
operator Pjv, which projects states onto a subspace of the 
Hilbert space with fixed total particle number N. The 
operator Pj^ itself does not preserve the normalization of 
a state and thus requires an additional rcnormalization 
factor Af. With this operator, we can show that any 
state of the form ([T]) can be written as a projection of 
the coherent product state 



M 



AI 



e 2 e'^'= vac) 



(11) 
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The number projection of this state reads 
PnUM) = e---^i' l^^l Vatc^" '3^^^ |vac) 

00 / M 

Ti=0 ' \ k 



N 



g 2 



I vac) 



(12) 



This coincides with the general inital many-body ground 
state (j?! with the normalization factor A/" = N'^ /Nl, 
which has the form of a Poisson distribution in the total 
particle number N. Note that we used |/3fcP = N. 

To perform calculations conveniently we can write the 
projection operator in the form of a phase integral 



1 



N 



J- 



dcf) e 



i{N-N)(l> 



(13) 



with the total particle number operator TV = X^i^l^^- 
Note that the projection operator P/v has the properties 
[Pn, bib,] = [Pjv, H] = for all z, j, and P^ = Pn/VM- 
The benefit of writing the BEC ground state as num- 
ber projection of the |{/3fe}) state, with the projector from 
Eq. ( 13 ) is that the SPDM then factorizes with additional 
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phase integrals that can be solved exactly. We consider 
the time-evolved SPDM under the evolution of Hamilto- 
nian Q, which is given by 



1 



e-'^*|BECjv) 

-iN(p 



2ttN . 

X {{PuW"'b\h,e-'"'e^-^\{Pk}) 



1 

27r7V 



-iNct> 



In the last line we have used the fact that 



(14) 



and 

In total we find (i ^ j): 



(21) 



(JV-l) 



71=0 



for every site k. 

We can then exactly evaluate the SPDM elements from 
Eq. ( 14 ) for time-evolution under the Hamiltonian Q . 
In the diagonal case, [H, blbi] = 0, and thus the diagonal 
elements are time-invariant. At t — 0, we find 



>j)t=0 



2TTj\f 



(16) 



where we used 



k/k 



l3k\Pk)k, and {a* 



\a\-/2-\i3\ /2_ Integrals of the form of Eq. (16) can 



be solved with the parametrization z = e^'^. This yields 



dz-^e 



Nz 



Af 



-, rr lim 

{N - 1)! 2^0 



d 



(17) 



where the integration takes place along the complex unit 
circle, denoted by C, and we again used |/3fcP = N. 
Thus, the initial SPDM factorizes into I3*l3j, and the di- 
agonal elements remain constant in time {b\bi)t = 

We finally calculate the time-evolution of the off- 
diagonal SPDM elements. We write the Hamiltonian Q 
as H = and find for i ^ j 



{{Pu}W"%b,e-^"'me^^}) 
X n {Pk\Pke'^)k. 



■iiiAt 



bje 



-ih^t 



(18) 



The factorized elements can be straightforwardly cal- 
culated as 

(19) 



where we again used the parametrization z — e'*^ and an 
integration along a complex unit circle. 



C. Revivals in a homogeneous system (periodic 
boundary conditions) 

In the case of a homogeneous density, /3i = f3 = 
\J N/M for all i. This implies that we assume periodic 
boundary conditions, and no external trapping (e^ = 
for all i). The off-diagonal SPDM elements then become 
site-independent, with Eq. (21 1 simplifying to 



N 



(2cos(C/t) - 2) 



(22) 



For N — M — 2, we reproduce the exact result from 



Eq. (10). In the limit where the particle number is much 

P, we find 



larger than the onsite density, i.e., N ^ 



N 



(2cos([/t) - 2) 



[N-l] 



|2g|/3|^[2cos((7t)-2] 



(23) 



since lim7v->oo(l + x/N) — . Hence, in this limit we 
reproduce exactly the result of the naive treatment from 
Eq. (If. 

For homogeneous densities, the height of the q = 
peak in the quasi-momentum distribution is given by 

n,.o = il/M)Ef!Mbj)t = {M- l){b%) + Thus, 
the experimentally measured visibility of timc-of-flight 
interference fringes is directly connected to the site- 
independent off-diagonal SPDM elements. In Fig. [T] we 
show a plot of the time-evolution of the these elements 
beginning from an ideal non-interacting ground state 
with density |/3p — 1, and including the first collapse 
and revival. For increasing particle numbers, we see that 
the results from Eq. (22 1 converge rapidly to the values 
given by the naive treatment in Eq. ([3|. Indeed, already 
for moderate numbers, N ^ 5, the exact results are very 
close to the approximate result. 
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FIG. 1. Plot of the time-evolution of the off-diagonal SPDM 
elements for an initial ideal superfluid state {Uo = 0) with 
homogeneous density and without external trapping [si = 0) 
in a deep optical lattice (J = 0). The density is chosen as one 
particle per site (|/3|'^ — 1), and we show the first collapse and 
revival for increasing number of particles. We compare results 
from the exact treatment of Eq. \22\ to the naive result from 
Eq. ([sjl (solid black line). Periodic boundary conditions are 
assumed. 
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FIG. 2. The time-evolution of the relative error et of the naive 
approach compared with exact results in the height of the 
quasi-momentum distribution peak (i.e. nq=o). The results 
are for a ID system with box boundary conditions, resulting 
in a site-dependent particle density. Results are shown for the 
time-evolution of the error until the first revival at t = 2n/U, 
and for increasing system sizes (A^ — M, Uo ~ 0, J — 0, 
= 0). 



We can define the relative error in the visibility of in- 
terference fringes, i.e. the relative error in the height of 
the nq=o peak between the naive and the exact treatment 
as 



'9=0 



'9=0 



(24) 



where n°^§*^ is obtained from Eq. (22) and n'^^?i'o° from 
Eq. (|3|. Using the parametrization at = |/3p(2 cos(J7i) — 
2), i.e. — 4|/3p < at < 0, and expanding in powers of 
1/7V, we find 



27V 



24Ar2 



(25) 



Thus, for large N and small homogeneous densities, the 
relative error from the naive treatment decreases propor- 
tional to Note that this error can also be written 
in terms of the system size M. Up to first order wc find 
for small densities 

2cos{Ut)-2 |/3|2(2cos(C/t) - 2)2 ,^ / 1 \ 
M 2M ^^[m^J- 

The error decreases proportional to 1/M. In the next 
section we will also treat the more general case of inho- 
mogeneous densities below. 



D. Revivals for inhomogeneous systems 

Inhomogeneous densities arise from external trapping 
potentials, i.e. from site varying energy-offsets Si in 
Eq. ^, or fr om confinement represented by open bound- 
ary conditions. The initial state for non-interacting par- 
ticles on the lattice can be found by solving the single 



particle problem for the amplitudes Pi. Here we will fo- 
cus on a ID lattice array with M sites, for which the 
initial single-particle Hamiltonian reads 



{\n){n + l\ 



en\n){n\, (26) 



where \n) denotes the state for a particle being at site 
n. Without external trapping (e„ = for all n) and 
for periodic boundary conditions {\M -I- 1) = |1)), the 
eigenstates of Hamiltonian ( 26 1 are simply given by 



IVm) — '^'jLi ^^P{'^1m.o-j)\j) / V M , with quasi-momenta 
Qm — 2iTm/aM and eigenenergies Em — — 2cos(gma). 

For the case of box boundary conditions, we can 
write the system on a grid of AI + 2 sites (la- 
beled by z = 0, ...,M + 1), and local energy offsets 
eo,eM+i -> 00, and e„^o,A/+2 = 0. To fulfill the 
boundary conditions, the eigenstates are then given by 
Wm) — X^^lo^ with normalization 

M = X^jio^ sin^(gmaj)- The quasi-momenta are given 
by Qm = rmr/{M + l)a, where m > 0, and the eigenen- 
ergies are E^. — — 2cos(gma). Thus, the coefficients f3i 
for the renormalized single-particle ground-state ([6| are 
given by 



ft 



M+1 



(27) 



In Fig. [2] we plot the relative error et between the 
exact and the naive treatment when we begin from an 
inhomogeneous initial state due to box boundary condi- 
tions. We show the time-evolution of the error in the 
height of the q = quasi-momentum peak, as defined in 
Eq. ( 24 ) , covering the time-scale of the first collapse and 
revival until t = 27r/C/ in a system with N = M parti- 
cles. We find that et generally decreases with increasing 
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FIG. 3. The scaling of the relative error et=Tr/4u of the naive 
approach compared with the exact results at the middle of the 
first collapse t = n/MJ . Results are shown both for periodic 
and open (box) boundary conditions. To compare the results 
to Eq. pSl, we also plot l/N {Uo = 0, J = 0, £i = 0). 



system size, and assumes its largest value approximately 
in the middle of the collapse process at t sa tt/AU. 

To analyze the scaling of the error in the particle num- 
ber N, in Fig. [3] we plot the relative error in the height 
of the q = quasi-momentum peak as a function of N. 
For homogeneous densities in the limit of large N, et is 
given by Eq. (25), i.e. it decays proportional to 1/iV. 

We thus find that for large N, the error scales as 
for both homogeneous and inhomogeneous densities stud- 
ied here. For typical ID lattice sizes in experiments of 
approximately 50 sites, the relative error becomes of the 
order of 1%. In this sense, the naive treatment provides a 
good semi-quantitative description of the system. How- 
ever, because experimental system sizes are of the order 
of 10-100 lattice sites, it should be possible to measure 
the quantitative discrepancies arising in the revivals as a 
result of finite system sizes and fixed particle numbers. 



E. Effects of finite tunnelling or initial interactions 

In this section we compare our results to numerical 
calculations for the ID Bose-Hubbard model, in order to 
account for the effects of a non-zero tunneling amplitude 
J during the time-evolution, and a non-zero initial inter- 
action strength Uq. We consider the time-evolution of the 
quasi-momentum-distribution in systems up to M = 51 
sites, using the TEBD algorithm [19l [20]. This algo- 
rithm makes possible the near-exact integration of the 
Schrodinger equation for ID lattice and spin Hamiltoni- 
ans based on a matrix product state ansatz. We optimize 
the algorithm for fixed total particle numbers as is done 
in Density Matrix Renormalization Group methods [26- 

In Fig. |4] we show results for the time evolution of the 
quasi-momentum distribution. We start with the exact 
ground state for box boundary conditions. It is obtained 
by an imaginary TEBD time-evolution and yields an ini- 




FIG. 4. Time-evolution of the quasi-momentum distribu- 
tion Uq for several tunneling amplitudes J in a system with 
A4 — 51 sites and box boundary conditions. Initially the 
system is in the exact ground state for Uq/Jq = 0. Pan- 
els (a)-(c) show exact results obtained numerically by using 
the TEBD algorithm. Panel (d) shows the evolution of the 
Uq^o peak. The TEBD results for J 



are compared to 
the exact results obtained from Eq. (21 1, and to the naive 
treatment obtained from Eq. ([2|. The particles are intitially 
non-interacting {Uo — 0). 



tial density profile, (3i, consistent with Eq. (27). The 



exact time-evolution calculation covers the first two re- 
vivals within the time < tU < 57r « 15.7. In panels 
(a)-(c) we show TEBD results for J = 0, J = OmU, 
and J = O.IU, respectively. In panel (d) we show the 
time-evolution of the peak of the quasi-momentum dis- 
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FIG. 5. Panel (a) compares the initial quasi-momentum dis- 
tribution for non-interacting particles {Uo ~ 0) with the one 
for weakly interacting particles {Uo/ Jo = 0.1). Shown are re- 
sults for a system with 51 sites and box boundary conditions. 
The inset shows the scaling of the n^^o peak as a function of 
Uo- Panel (b) shows a comparison of the time-evolution of the 
Uq^o peak for initially non interacting particles (solid-lines) 
with weakly interactings ones {Uo/ Jo = 0.1). 



tribution, i.e. Uq^Q. There, we compare the results from 
our TEBD calculations with the exact calculations for 
0, i.e. from Eq. (21), and the naive treatment from 



J= 

Eq. ([2]). For our numerical calculation, convergence in 
several numerical parameters had to be checked. For the 
results obtained in this section we have tested conver- 
gence of the TEBD algorithm for numerical parameters, 
and use matrix sizes up to x = 200 [HI [50] , and allow 
occupation numbers up to 14 particles per lattice site. 

For J = results from the analytical number projec- 
tion method (Eq. (21 1) perfectly coincide with the TEBD 
results as expected, whereas visible differences are ob- 
served to the results from the naive treatment. These 
errors are small and scale as in Fig.[3| For non-zero tun- 
neling, we find that the revivals are strongly damped, 
i.e. the peak in the momentum distribution for g = is 
smaller, and the quasi-momcntum distribution at the re- 
vival time is broadened compared to the initial one. In 
addition, we see that the point in time of the first and 
the second revival are markedly shifted towards earlier 
times for increasing tunneling strength J. Both effects 
that appear here in exact numerical calculations of an 
ID systems have been recently calculated perturbatively 



|21) . and for hard-core bosons as well as on a mean-field 
level in high dimensions [22j . 

For non-zero initial interaction strength Uq > 0, we 
note that the initial quasi-momentum distribution is 
slightly modified. However, the basic structure of the re- 
vivals is unaffected for Uq Jo. This is shown in Fig. [s] 
for a system of Af = 51 sites and box-boundary condi- 
tions. We find that due to boundary effects for Uo < 
O.lJo, the quasi-momentum distribution for finite Uq be- 
comes narrowed compared to the one for non-interacting 
particles (Fig. [sf^a)). Further increasing the initial inter- 
action strength leads to a depletion of the q = state 
and a broadening of the quasi-momentum distribution 
(see inset). In Fig. [sjb) we compare the time-evolution 
for initially non-interacting particles (solid lines) with the 
one for Uq/Jq = 0.1. We show Uq^Q for tunneling ampli- 
tudes J/U — 0, 0.05, and 0.1 during time-evolution. We 
find that the point in time of the first and the second re- 
vival is not altered significantly due to C/q 7^ 0. However, 
we find that the structure of the peaks modifies and is 
smoothed for J/U > 0. 



IV. SUMMARY & OUTLOOK 

By using a number projection method, we have shown 
that exact results can be derived for the collapse and re- 
vival of the SPDM for atoms in an optical latittce. We 
properly account for a fixed total number of particles in 
the system, and are able to take into account the effects 
of external potentials. For large numbers of particles N, 
we find that the discrepancy to naive calculations, based 
on assuming a coherent state on each site in the particle 
number scales as 1/iV. For as few as 5 particles on 5 sites, 
this approach can provide a good semi-quantitative ap- 
proximation. However, the discrepancies should be mea- 
sureable in experiments. Effects of finite tunnelling in 
the lattice, which can be treated perturbatively [21], or 
exactly for hard-core bosons [22] can be investigated di- 
rectly via time-dependent numerical calculations. 

The general technique we use here could be extended 
to the case in which the onsite interaction strength is de- 
pendent on the occupation number, as has been observed 
in recent experiments ^Tj. 
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